Optical and dc conductivities of cuprates: 
Spin-fluctuation scattering in the t-J model 
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A microscopic theory of the electrical conductivity o~(ui) within the t-J model is developed. An 
exact representation for <r(w) is obtained using the memory-function technique for the relaxation 
function in terms of the Hubbard operators, and the generalized Drude law is derived. The relaxation 
rate due to the decay of charge excitations into particle-hole pairs assisted by antiferromagnetic spin 
fluctuations is calculated in the mode-coupling approximation. Using results for the spectral function 
of spin excitations calculated previously, the relaxation rate and the optical and dc conductivities 
are calculated in a broad region of doping and temperatures. The reasonable agreement of the 
theory with experimental data for cuprates proves the important role of spin-fluctuation scattering 
in the charge dynamics. 

PACS numbers: 72.10.-d 72.10.Bg 78.20.Bh 71.27. +a 



I. INTRODUCTION 

Studies of charge dynamics in superconducting 
cuprates provide valuable information concerning elec- 
tron interaction with bosonic modes which is impor- 
tant for elucidating the pairing mechanism in high- 
temperature superconductors. There is a vast literature 
devoted to these studies (for reviews see, e.g., [lJ-Q) ■ 
Two major scenarios have been proposed, the electron- 
phonon coupling and electron interaction with the an- 
tiferromagnetic (AF) spin fluctuations. Angle-resolved 
photoemission spectroscopy (ARPES) points to an im- 
portant role of spin fluctuations in a renormalization of 
the single-electron excitation spectrum (see, e.g., Q and 
references therein) which is supported by measurements 
of the infrared (IR) absorption in a wide region of temper- 
atures and doping (see, e.g., 6-8] and references therein). 
The main argument against the spin-fluctuation pair- 
ing mechanism, a weak intensity of spin fluctuations at 
the optimal doping seen in inelastic magnetic neutron 
scattering experiments [9(, was dismissed in recent res- 
onant inelastic x-ray scattering [lCJ. In a large family 
of cuprates AF paramagnon excitations with dispersions 
and spectral weights similar to those of magnons in un- 
doped cuprates were found. However, a decisive role of 
the electron-phonon interaction (EPI) has been claimed 
in a number of theoretical studies (for a review see [HI). 

The optical conductivity (OC) <t(lo) of cuprates re- 
veals a complicated evolution with doping and tempera- 
ture. The undoped parent compounds are AF insulators, 
where the OC exhibits a peak at the charge-transfer en- 
ergy lo < 1.8 eV. Under doping, the insulator-to- metal 
transition occurs when the charge-transfer gap is filled 
up with states and the spectral weight is transferred to 
the lower energy, the Drude peak at lo — > with a width 
to < 600 cm -1 and a broad mid-infrared (MIR) band 
at higher energies, lo < 5000 cm -1 . With increasing 
hole concentration, the MIR absorption shifts to lower 



energies and merges with the Drude contribution (see, 
e.g., Refs. [12H16I ]). The Drude peak significantly nar- 
rows with decreasing temperature and is attributed to 
the relaxation of coherent quasiparticles, while the origin 
of the temperature independent MIR contribution is still 
under discussion. In Ref. [17| the metal-to-insulator tran- 
sition (MIT) was studied by measuring the OC for the 
underdoped Bi-based and YBCO-based compounds for 
hole concentrations from 6 = 0.12 to S = 0.015. With 
decreasing hole concentration, the Drude peak at low 
temperatures transforms into a far-infrared (FIR) band 
at energies lo > 200 cm -1 which acquires a gap at the 
MIT for hole doping S < 0.07. Note that the onset of 
the metallic phase occurs at a doping much higher than 
S ~ 0.02 at which the AF long range order (LRO) van- 
ishes [Hj]. 

It is generally believed that superconducting cuprates 
are doped Mott-Hubbard (charge-transfer [l9|]) insula- 
tors, where the insulating phase of the undoped parent 
compounds appears due to a strong Coulomb repulsion, 
Hubbard U > U C 2, where the critical value U C 2 for the 
MIT is larger than the electronic bandwidth W (see, 
e.g., [13]). In this case the AF LRO in the undoped com- 
pounds originates from the strong AF exchange interac- 
tion characteristic for Hubbard systems. However, it is 
also possible to explain the insulating phase as caused by 
the AF energy gap induced by the AF LRO where the 
Coulomb interaction plays a secondary role. In recent 
publications, this problem was discussed by analyzing the 
OC for typical electron-doped Nd 2 _ 2; Ce a; Cu04 (NCCO) 
and hole-doped La2- a; Sr a ;Cu04 (LSCO) compounds. In 
Ref. [2JI, the OC was calculated for the Hubbard model 
in the paramagnetic and the AF phases. Using the dy- 
namical mean-field theory (DMFT) [22j , the optical spec- 
tral weight (given by the restricted sum rule with the 
integration over energy in Eq. ^ up to fl = 0.8 eV) 
was calculated. Comparing the doping dependence of the 
theoretical and measured spectral weights for the NCCO, 
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LSCO and other cuprate compounds, it was found that, 
in the paramagnetic phase, the fitted U is smaller than 
the critical value U C 2 ~ 1.5 W for yielding the Mott- 
Hubbard insulating phase. At the same time, the AF 
phase provides the insulating state for the undoped sys- 
tem at the fitted U. So, it was concluded that antiferro- 
magnetism is essential in producing the insulating state. 
A different conclusion concerning the hole-doped 



cuprates, was obtained in Refs. [23|, l24j . The OC was 
calculated for a realistic three-band p-d model for NCCO 
and LSCO using the local density approximation com- 
bined with the DMFT [25]. It was found that, whereas 
for the electron-doped NCCO the AF interaction is neces- 
sary to yield the insulating undoped state, the hole-doped 
LSCO belongs to the Mott-Hubbard system, where the 
insulating state is due to strong electron correlations but 
not to the AF interaction. However, in the DMFT the 
short-range AF correlations are neglected. As shown in 
Ref. 26], by taking them into account, we find a much 
narrower Hubbard band W which leads to a lower criti- 
cal parameter U C 2 that may put the electron-doped com- 
pounds also into the class of Mott-Hubbard insulators. 
Therefore, for a more accurate estimation of U C 2 the clus- 
ter DMFT [13 should be used. By this method the OC 
was calculated in the t-J model in Ref. [28|. The results 
reproduced quite well the delicate changes of the spectral 
weight transfer at the transition to the superconducting 
phase observed in experiments (see, e.g. [29|, [3(| and ref- 
erences therein). This proves that the t-J model cap- 
tures the essential physics of the low-energy excitations 
in cuprates. 

In the limit of strong correlations, extensive numeri- 
cal studies of the OC within the Hubbard model and t-J 
model for finite systems have been carried out. Earlier 
results are reported in Refs. (3ll - [35l | . Due to a small size 
of the clusters, only one- or two-hole motion was consid- 
ered. In that case the Drude peak cannot be resolved, 
and its intensity versus doping was studied by calcula- 
tions of the kinetic energy (see Eq. (f20|) L Several peaks 
found in the OC in the MIR region may be related to 
local spin excitations at u> ~ 2 J and to string excitations 
at higher energies. This observation was confirmed in 
analytical studies of the charge correlations in a weakly 
doped t-J model using the cumulant expansion within 
the Zwanzig-Mori projection technique [361 ] . The peaks 
found in the OC in the energy region u) ~ 2 J were as- 
signed to excitations due to internal degrees of freedom 
of the spin-bag quasiparticles. 

In several studies, an important role of the EPI result- 
ing in polaronic effects was stressed. However, contra- 
dicting explanations were proposed for the two absorp- 
tion bands observed in the MIR region, one at a lower 
energy near the FIR region and another at higher en- 
ergies. In Ref. |37j . the OC of one hole in the Holstein 
t-J model using the DMFT was calculated. It was shown 
that the IR absorption is characterized by the coexistence 
of a magnon peak at low energy and a broad polaronic 
band at higher energy. The two absorption bands were 



explained in Ref. [38[ by the coupling of a hole to two 
kinds of bosonic excitations. The lower energy peak at 
qj ~ 1000 cm -1 was ascribed to the phonon sideband, 
while the higher energy peak at ui ~ 4600 cm -1 was 
considered as the magnon sideband of the lower peak. 
As discussed in Ref. [39(, the two-peak structure in the 
MIR region may be explained by the coupling of a doped 
hole to magnetic excitations. The low-energy peak repre- 
sents the local magnetic excitation, attached to the hole, 
while the higher-frequency peak corresponds to the MIR 
band that originates from coupling to spin-wave exci- 
tations, broadened and renormalizcd by phonon excita- 
tions. Thus, the studies of the Holstein t-J model suggest 
that the complicated absorption structure experimentally 
found in the MIR region is caused by magnetic excita- 
tions which are coupled to phonons via doped holes. 

In analytical studies, the OC is frequently evaluated 
within the simple electron-hole (bubble) diagram approx- 
imation for the current-current correlation function pro- 
posed by Allen [40l ]. The finite-temperature version of 
the Allen approximation was derived in Ref. j4l| . This 
method was used in studies of the optical IR data within 
electron-phonon models (see [ll| and references therein) 
and spin-fermion models (see, e.g., [42|,|43[ and references 
therein). 

A general approach based on the Mori memory- 
function method 441 for the calculation of the current- 



current relaxation function was proposed in Refs. 45[. In 
this method, the transport relaxation rate is expressed 
directly in terms of the force-force relaxation function 
which can be further evaluated by perturbation theory 
with a proper consideration of the wave-vector depen- 
dence of the transport vertex. In particular, in Ref. (46| 
the memory-function method was used to calculate the 
OC in the limit of strong electron correlations within the 
Emery model for Cu02 plane. The relaxation rate for 
electrons scattered by AF spin fluctuations was calcu- 
lated in a fair agreement with experiments. In Ref. [47j |. 
by taking into account only the incoherent part of the 
electronic spectrum, a scaling expression for the fre- 
quency dependence of the relaxation rate and the con- 
ductivity in the t-J model was obtained. 

Electron interactions with bosonic modes can be re- 
vealed in the low-energy part of the OC and dc conduc- 
tivity. To shed more light on the scattering mechanism 
in cuprates, in the present paper we calculate the op- 
tical and dc conductivities within the t—J model. The 
main goal of this work is to demonstrate that AF spin- 
fluctuation scattering is the essential mechanism of the 
low-energy charge dynamics of underdoped and optimally 
doped cuprates. 

Using the memory-function method, we derive an 
equation for the relaxation rate which is determined by 
the kinematic interaction for the Hubbard operators and 
depends only on the parameters of the t-J model: the 
hopping parameters and the AF exchange coupling. The 
relaxation rate is calculated by taking into account elec- 
tron scattering by spin fluctuations which are described 
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by the spin-excitation spectral function calculated in our 
previous works [H, H^] . Therefore, we are able to con- 
sider effects of spin excitations on the charge dynamics 
within a microscopic theory without fitting parameters. 
In our approach we obtain a reasonable agreement with 
experiments for the relaxation rate, the optical conduc- 
tivity, and the resistivity in broad regions of temperature 
and doping, in particular, in the underdoped region with 
a strong AF short-range order (SRO). 

In Sec. II we formulate a general theory of the opti- 
cal conductivity within the memory-function formalism. 
The application of this theory to the t-J model is given in 
Sec. III. Numerical results and discussion are presented 
in Sec. IV. In Sec. V we summarize our results. 



II. MEMORY-FUNCTION THEORY 

In the linear response theory of Kubo [50], the fre- 
quency dependent conductivity is defined by the current- 
current relaxation function, 



Oxx(u) = — ((J x \J x )) a 



1 

VJ 



dte iuit (J x (t),J x ), (1) 



where V is the volume of the system. Here, the Kubo- 
Mori scalar product 



(A(t),B) = / d\(A(t - i\)B) 



(2) 



for the operators in the Heisenberg representation, 
A(t) = exp(iHt)Aexp(—iHt), is introduced . (AB) de- 
notes the equilibrium statistical average for a system with 
the Hamiltonian H (here f3 = l/T, h = k B = 1). The 
real part of the conductivity ([IJ obeys the sum rule 



7T 17T 

duRea xx (u) = —Xo= 2y([Jx,Px})- (3) 



Here P x = eJ2i^i^i is t ne polarization operator. Rf 
is the x-component of the lattice vector pointing to site 
i, e is the electron charge, and Ni is the number opera- 
tor. The current operator is defined by the time deriva- 
tive of the polarization operator: J x (t) — dP x (t)/dt = 
P x (t) = —i[P x ,H]. The static current-current suscepti- 
bility xo = (Jx, Jx) is related to the effective number of 
charge carriers participating in the absorption, 



Neff = 



2mvo 



ire 



diu~R,ea xx (u) 



2 N 



Xo ■ 



(4) 



where vq = V/N is the unit cell volume and m is the free 
electron mass. The sum rule Q is frequently written 
in terms of the plasma frequency w p i defined by: ui^ — 
^Xo/V = U)l M N cS , where w£ pl = 4nNe 2 /mV. 

To calculate the conductivity ([T]), it is convenient to 
employ the memory-function approach of Mori [44j by 



introducing the memory function M(lo) for the relaxation 
function 451 . 



= (J x \ Jx)\ 



Xo 



M(w) 



(5) 



From the equations of motion for the relaxation func- 
tion $(i — t') = ((J x (t)\J x (t')j) the memory function is 
determined by (see Appendix A) 



M(u) = ((F x \F x ))t° pel) (l/Xo), 



(6) 



where F x = iJ x = [J x , H] is the force operator. The defi- 
nition of the memory function ((6]) as the "proper" part of 
the force-force relaxation function is equivalent to the in- 
troduction of the projected Liouvillian superoperator for 
the memory function in the original Mori technique [Hf. 

Using Eq. ([5]), the frequency dependent conductivity 
(JI]) can be written in the form of the generalized Drude 
law: 



crxx(w) = a(ui) 



An 



m 



1 



m(uj) T{u)-iu' 



(7) 



where the effective optical mass and the relaxation rate 
are given by 



m(uj) 
m 

AM 



= i + a(w), r(w) = 

= M'(u))/u, r(w) = 



1 + A(w) : 
M"{uj). 



(8) 
(9) 



Here the real and imaginary parts of the retarded mem- 
ory function M(lu + i0 + ) = M'(cj) + iM"(u>) are intro- 
duced . They are coupled by the dispersion relation: 



M'(w) = - 



00 , M>') 
dcu — ■ . 



(10) 



Both the real and imaginary parts of the memory func- 
tion can be directly related to experimental data for the 
inverse conductivity ([7]) 0: 



r( w ) = ^lRe— -. 
47T cr{oj) 



i + AM = - r eLim- rT . 

47tu; cr(wj 



(11) 



Below we calculate these functions for the t-J model. 



III. RELAXATION RATE 



We consider the t-J model on the square lattice which 
in the conventional notation reads: [U [52| 

H = H t + Hj = -J2 Ufaj a ~a ]a - fij^ N i 



(12) 



where ty is the hopping integral and is the AF ex- 
change interaction. Here af^ = af a (l — n is ) is the pro- 
jected electron operator with spin a/2 — ±1/2, (a = 
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—a) on the lattice site i, Ni = J^a 'ha^ lB tne num- 
ber operator, and Sf — (1/2) J2 S s > ®ts a s a '°»«' ^ s ^ ne a " 
component of the spin operator (erf s , are Pauli matrices). 
The chemical potential \i is determined from the equation 
for the average electron occupation number (Ni) = 1 — 5, 
where 5 is the hole concentration. 

To take into account the projected character of elec- 
tron operators we employ the Hubbard operator (HO) 
technique [H H]. The HO xf = \i,a)(i,/3\ de- 
scribes the transition from the state \i,/3) to the state 
\i,a) at the site i, where a and /3 denote three pos- 
sible states: an empty state (a, /3 = 0) and a singly 
occupied state (a, f3 = a). The completeness relation 
X ( 00 + J^cr Xf a = 1 rigorously preserves the constraint 
of no double-occupancy of any lattice site. From the 
multiplication rule xf = 5 M Xf follow the com- 
mutation relations for the HOs: 



[Xf ,X] s ] ± =5 t3 ( 



s P-yXf s ± 1 



10 



(13) 



The upper sign refers to the Fermi-like operators creating 
(Xf°) or annihilating (X® a ) electrons, while the lower 
sign refers to the Bose-like operators, such as the number 
or spin operators: 



Ni — ^ X™, S? — - a Xf CT , Sf — X° 



(14) 



The commutation relations result in the kinematic inter- 
action for HOs (see Eq. (|B3| ). Note that the term "kine- 
matic interaction" was introduced by Dyson (55j for spin 
operators. 

Using the HO representation for af a = Xf° , a Jcr = 
X? a , and Eq. CE 
form: 



we write the Hamiltonian (IT21 in the 



H = - y, ti 3 x?xy-nY, x t 



- Y JiAxrxr-xrx 



era \ 
3 ) 



(15) 



The relaxation rate T(uj) = M"(uj) is calculated by 
Eq. ^ in the mode-coupling approximation (MCA) as 
described in Appendix B. In this approximation we ob- 
tain: 



r( w ) 



(e^-l) 27re 2 



N 



E 



duJiduj2duj^ 



xn(wi)[l — n(uj2)]N (ws) S(u>2 — usi — ^3 + oS) 
x^(k,k-q)x^(q,W3)A(k,a;i)A(k-q,W2) , (16) 

where n(u>) = (exp/Sw+l) -1 and N(uj) = (exp/3cj— l) -1 . 
The momentum dependent (transport) vertex is given by 

g x (k,k-q) = v x (k) t(k - q) - ^(k - q) t(k) 

- J(q)/2K(k)-« x (k-q)], (17) 



where t(k) and J(q) are the Fourier transforms of 
the hopping integral and the exchange interaction, and 
v x (k) = —dt(k)/dk x is the electron velocity (see Ap- 
pendix B). The spectral function for the charge-spin ex- 
citations XcsiQi w ) = (l/ 7r )I m Xcs(q, u) is defined by the 
corresponding commutator Green functions (GFs), 



Xcs (q,c) = -(l/4)«iV q |iV_ q )) u - «S q |S_ q )) U) 



(18) 



where we used Zubarev's notation [56J for the retarded 
two-time GFs. The spectral function of electronic exci- 
tations is defined by the imaginary part of the anticom- 
mutator electronic GF, 



A(k, W ) = -(l/7r) lm((X°°\X^ 



(19) 



The static current-current susceptibility xo is connected 
with the effective number of charge carriers Q which for 
the t-J model reads 



m 



N «* = '-fYK R i- R ^ 2t ^ x i Ux J a ) 

m x ^ d 2 t(k) /vaQ v0ax 
-N^^kT' ( k k } " 

k,cr x 



(20) 



For the t-J model with the nearest-neighbor hopping 
only, (Rf — R X j) 2 = a 2 , where a is the lattice param- 
eter, the effective number of carriers is related to the 
averaged kinetic energy, N e g = (a 2 m/N)(—H t ) . This 
relation is often used in the calculation of the charge 
stiffness (Drude weight) in finite-cluster studies (see, e.g., 
Refs. (HUGH). 



IV. RESULTS AND DISCUSSION 

A. Spectral functions 

In numerical calculations we have to use models for the 
charge-spin susceptibility (TTSJ) and the one-electron spec- 
tral function (1191) . The spin-excitation contribution in 
Eq. (TT8")) is described by the spectral function x"(q,w) = 
(3/2)x£(q,w), where x£(q,w) = -(V^Im <<S+|SlJL. 
For the latter we use the function calculated in Ref. [48[ 
for the t-J model, 



x±(q,«) 



■wV(q) (m(q)/7r) 



■wS,'(q,w) 



>5V'(q)] 



.(21) 



Here, the spectrum of spin excitations in the generalized 
mean-field approximation w q determines the static spin 
susceptibility x q = m(q)/u; q with m(q) = ([iS q , 5T q ]), 

where iS+ = [S+,H] . The self-energy E s (q, w) = 
E s '(q,w) + iE s "(q,w), where E/(q,w) and E s "(q,w) 
are the real and the imaginary parts, respectively, is 
determined by the many-particle relaxation function 
E s (q,w) = (l/m(q)W-5+| - #I q ))P roper calculated in 
MCA (see Refs. [48|, |49(). Taking into account that the 
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main contribution to the relaxation rate (116)) from the 
spectral function (|2"Tj) is given by frequencies close to 
the renormalized spin-excitation frequency w q = [w^ + 
w q E^,(q, Wq)] 1 / 2 , we approximate the damping of spin 
excitations by the function S s "(q) = E/'(q, ui = u; q ). 

To calculate the contribution to the relaxation rate 
(fT6|) from charge (density) excitations in Eq. (fT8|) . we use 
the spectral function calculated in Ref. [57j for the t-J 
model. Our results show that charge excitations give the 
main contribution in the region of high energies, uj ~ it, 
which, however, is several times weaker than the spin- 
excitation contribution and, therefore, can be safely ig- 
nored. The different energy scales for spin excitations, 
lj ~ J, and density excitations, u> ~ t, were also found in 
an exact diagonalization study of the t-J model [34j . 

The self-consistent solution of the system of equations 
for the spectral function (fT9|) and the single-electron self- 
energy in Ref. [HI has shown that close to the Fermi en- 
ergy there appear well-defined quasiparticle excitations. 
This result permits to approximate the spectral function 
by the expression 



A(k,u) =Q5(u}-e k ), 



(22) 



where Q = 1 — rt/2 is the spectral weight for electronic ex- 
citations in the t-J model. To model a realistic electronic 
spectrum which crosses the AF Brillouin zone (BZ), as 
observed in ARPES experiments (see, e.g., Ref. [531), we 
consider the model dispersion 



-4Q[t ai 7 (k) + t'a 2 7' (k) + t"a 2 j"(k)} - p, (23) 



where t and t' — O.lt, t" — Q.2t are the hopping param- 
eters for the nearest and further-distant neighbors, re- 
spectively, and 7(k) = (l/2)(cosak x + cosak y ), 7'(k) = 
cos ak x cos ak y and 7"(k) = (l/2)(cos2afcj, +cos2afcj,). 
The kinematic interaction for the HOs results in a renor- 
malization of the spectrum (|23l) determined by the pa- 
rameters ai = [1 + Ci/Q 2 ] and a 2 = [1 + C 2 /Q 2 ], 



(SiS i±aa! / ay ) and C 2 



(SiSi 



where C\ 

(SiSi-j- 2a!C /2a B ) are the spin correlation functions for 
the nearest and the second neighbors, respectively (see 
Ref. |58j). With increasing doping, the effective band- 
width W of the dispersion (|2"3"|) increases due to the de- 
crease of AF SRO described by the doping dependence of 
the spin correlation functions in the renormalization pa- 
rameters a%, a 2 . In particular, for S = 0.09 (0.2) at T = 

we have W = 1.14 1 (2.78 1) in comparison with the un- 
renormalized bandwidth W = 8 1 Q = 4.36 1 (4.8 1). 

The Fermi surface (FS) determined by the equation 
£k F = is shown in Fig. Q] for various doping. The 
renormalization of the spectrum induced by the AF SRO 
provides a FS with hole pockets at low doping which is 
equivalent to a pseudogap in the (±tt, 0) and (0, ±7r) re- 
gions of the BZ. In the study of the electronic spectrum 
and the FS of cuprates, a more accurate calculation of 
the spectral function (|22l) including the self-energy con- 
tribution was performed as reported in Refs. (58j and [60j 
for the t-J model and in Ref. I2a| for the two-subband 



Hubbard model. In particular, in Ref. [26j the spectral 
function on the FS close to the (71", 7r)-point of the BZ has 
a weak intensity resulting in an arc-type FS as observed 
in ARPES experiments. By taking into account that the 
spectral function of spin excitations (|2"Tj) is peaked at the 
AF wave vector Q = (ir, it) and is very broad in other 
parts of the BZ (see Ref. [49j|). in the calculation of the 
relaxation rate (TTBl . only those parts of the FS are im- 
portant which are coupled by the AF wave vector Q. 
Therefore, the parts of the FS in Fig. [T] far away from 
the AF BZ, in particular the back side of the hole-pocket 
near (w/2,tv/2) of the BZ, give small contributions at 
the integration over the BZ in Eq. (TTBl . This reasoning 
justifies the quasiparticle approximation (|2"2"|) used in the 
calculation of the relaxation rate. 




(0, 0) 



FIG. 1: Fermi surface at 5 = 0.09 (solid line), 0.16 (dashed 
line), and 0.2 (dots). 



The effective number of charge carriers (j20j) is conve- 
nient to write in the form 

N oS (6) = t}K(S) = r)Q n(i k ) [cos(ak x ) (24) 



k.<r 



+ 2(t'/t) cos(ak x ) cos(ak y ) + 4(i" /t) cos(2ak x )]. 

For the spectral function (|2"2")l the average electron occu- 
pation number is given by (X£°X^ a ') = Q n(ek) ■ The 
prefactor 77 = 2ma 2 t = t (p 2 /2m)~ 1 is a dimensionless 
ratio of the hopping parameter t to the kinetic energy 
of an electron with the momentum p. In particular, for 
a = 3.8 A and t = 0.4 eV we have 7? = 3.79 t [eV] = 1.52. 

In these approximations, the relaxation rate (1161) is 
determined by the expression 



r(w) = t 



ir(e^ - 1) 3t 2 Q 2 



ujK(6) 2N 2 



/OO 
^Vi(qV) 

K,q 00 



x (^(k, k - q) <5(e k - e k - q +u)' -u) 
xJV(a/)n(ek)[l-n(E k _ q )], 



(25) 



where, using Eq. (|17p . the dimensionless transport vertex 
g x (k, k — q) = (1/at 2 ) ^(k, k — q) is introduced. 



G 



The real part of the conductivity (J7J may be written 

as 

^ =A ^ 5A FT^ra¥ - (26) 

where A = w 2 pl / (47rf) = e 2 / (muo*)- 

The real part of the memory function M'{uj) = uj A(w) 
is calculated by the dispersion relation (fTUj) using the 
relaxation rate (f25j). This enables us to calculate the 
effective optical mass ([8]), fh/m = 1 + A(0). 

In numerical calculations we take J = 0.3 1 and t = 
0.4 eV as an energy unit (0.4 eV = 3226 cm" 1 = 4640 K). 
The results for the relaxation rate and the optical conduc- 
tivity as functions of frequency, temperature, and hole 
doping are in a good overall agreement with experiments. 
This will be detailed in the following. 



B. Relaxation rate 

At zero frequency, the relaxation rate is related to 
the dimensionless electrical resistivity p = 1/ct(0) by 
r(0)/t = Neffp. The temperature dependence of p for 
S = 0.09, 0.16, and 0.20 is shown in Fig.dJ For a doping 
near and larger than the optimal doping (5 = 0.16), we 
obtain a nearly linear temperature dependence, as is also 
observed in experiments (see, e.g., Refs. [6ll. l62l|). 

The effective number of charge carriers iV c g given by 
Eq. is shown in the inset of Fig.[5J It does not reveal 
a notable temperature dependence and can be approxi- 
mated by the function N c s ~ 2 5. It is remarkable that 
N e g increases more rapidly than the hole doping. This 
result is in agreement with the in-plane optical conduc- 
tivity data on LSCO compounds (Refs. 0, EH) which 
yields the effective number of charge carriers N e g(uj) in- 
volved in optical excitations up to the cut-off frequency 
lu (upper limit of the integral in Eq. (|4])). N e g was found 
to be nearly proportional to 26 for doping S < 0.15, e.g., 
at 6 = 0.1, N cf f(u = 1.5 eV) = 0.19 (Fig. 11 in Ref. Qj) 
and N cS (u} = 2eV) = 0.26 (Fig. 10 in Ref. pj]). 




FIG. 2: Resistivity p(T) = 1/5(0, T) for doping 5 = 0.09 
(solid) , 0.16 (dashed), and 0.2 (dotted). In the inset the 
effective number of charge carriers N c h(S) at T — is shown. 



The frequency dependence of the relaxation rate T(lj) 
(|25[) at different temperatures and doping is plotted in 
Fig. [3] We obtain an increase of T(u>) with increas- 
ing temperature, which qualitatively agrees with exper- 
iments. In the overdoped case, the relaxation rate de- 
creases resulting from the suppression of spin fluctua- 
tions. The broad maximum in the frequency dependence 
of r(w), clearly revealed at low doping in Fig. [3] (a), 
shifts to higher frequencies with increasing doping. The 
doping-dependent hnite effective bandwidth W limits 
the highest frequency for the relaxation, ui < 2W, so 
that at very high frequencies, T(u>) vanishes according to 
r(cj) oc 1/uj — > 0. Let us point out that a maximum in 
the relaxation rate is also observed in experiments for the 
underdoped samples as, e.g.,_in YBa2Cu30 y (YBCO^) at 
cj - 2000 cm" 1 for y < 6.5 [13. 

The real part of the memory function M'(ui) shown in 
Fig. H] exhibits a maximum which height decreases with 
increasing temperature and doping. But the energy of 
the peak does not change with temperature as observed 
in experiments (see, e.g., [sHl])- In the underdoped case 
(S = 0.09), the temperature dependence of M'(ui) is very 
strong, as compared with the overdoped case (5 = 0.2). 
This results from the strong AF SRO at low doping that 
strongly depends on temperature. With increasing dop- 
ing, both the SRO and the influence of temperature are 
weakened. Qualitatively, the relaxation rate (Fig. [3]) re- 
veals the same trend. 

In Fig. [5] we show the temperature dependence of the 
effective optical mass at zero frequency, fh/m = 1 + A(0), 
at various doping. At small doping, a strong temperature 
dependence of fh/m is observed, which may be explained 
similarly as for M'{oj). For the overdoped case, the effec- 
tive mass shows a weak renormalization, rh/m ^2. In 
numerical studies of the one-hole motion in small clus- 
ters, a weak renormalization of the optical mass was de- 
duced at high temperatures, T > 0.2t (Ref. (35[) which 
is in agreement with our results. In early experiments, 
Ref. [12J, a large effective mass renormalization was ob- 
tained in LSCO ranging from rh/m ~ 25 for 5 = 0.1 
to fh/m = 16 for S = 0.2. However, later on, e.g., in 
Ref. [16| , a nearly doping independent modest renormal- 
ization of the effective mass, fh/m = 3 — 4, was observed 
in both LSCO and YBCO which is close to our finding. 

C. Optical conductivity 

The frequency dependence of the conductivity (|26[) 
for various temperatures and hole doping is shown in 
Fig. IH The temperature and doping dependences of 
the conductivity show a reasonable qualitative agreement 
with experiments (see, e.g., Ref. [III). At low frequen- 
cies, a large Drude peak is found which significantly 
narrows at low temperatures. We also obtain a broad 
and nearly temperature-independent MIR maximum at 
oj < 2t ~ 6000 cm.- 1 which slightly shifts to lower fre- 
quencies and becomes of lower intensity with increasing 




FIG. 3: Temperature dependence of the relaxation rate T(uj) 
for (a) S = 0.09, (b) 0.16, and (c) 0.2. Note the different 
scales. 



FIG. 4: Temperature dependence of the real part of the mem- 
ory function M'(u) for (a) S = 0.09, (b) 0.16, and (c) 0.2. 
Note the different scales. 



doping as observed in experiments (see, e.g., |16|). In our 
theory, the MIR absorption results from electron interac- 
tion with spin fluctuations which influence the electron 
scattering so that it decreases with increasing doping. 



D. Quantitative comparison to experiments 

Let us first compare the resistivity p(T) — (l/A)p, 
where A is given in Eq. ([26]), with experimental^ data for 
the underdoped cuprate YBCOg.5 (8 — 0.09) 6] shown 
in Fig. [7J (a). Here we use the value A = e 2 /(mvot) = 
8.00 • lO^ft-cm]- 1 taking v = 57.7 (A) 3 [6]. With- 
out a fitting procedure, we obtain a remarkably good 
agreement with experiment, both in the absolute values 
of the resistivity and in its temperature dependence. In 
Fig. [7] (b) we compare the resistivity p{T) with experi- 
mental data on LSCO [6l[ for the underdoped (8 = 0.08) 
and nearly optimally doped (8 — 0.17) samples with our 
results. The value of A = 4.86 • 10 3 [^-cm] _1 is ob- 
tained using vq — a 2 d — 95.3 (A) 3 (a — 3.8 A and 
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5 = 0.16 

5 = 0.2 
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FIG. 5: Temperature dependence of the effective optical mass 
fh/m = 1 + A(0) at various doping. 

d = 6.6 A). A reasonable agreement is observed at high 
temperatures, while at low temperatures our values are 
much smaller. An additional scattering mechanism, e.g., 
impurity scattering, should be invoked to explain the ex- 
perimental data. The comparison of our results with op- 
timally doped and overdoped Tl-compounds 63] shows 
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FIG. 6: Temperature dependence of the optical conductivity 
ff(w) for (a) S = 0.09, (b) 0.16, and (c) 0.2. 



the same trend. 



Now we estimate the plasma frequency w p i = 

[Ncg] 1 ^ 2 wo, pi- For the optimally doped case, 8 = 0.16, 
we have N e g — 0.3 and u p \ = 0.55wo, P i ■ For LSCO 
( WO pl = 3.72 eV) we get uj p i = 2.05 eV and for YBCO 
(w ,pi = 3.96 eV), lu p \ = 2.18 eV. These values are close 
to experiments, while the LDA calculations in Ref. [13] 
give the somewhat larger value ui p \ ss 2.9 eV. For the un- 
derdoped YBCO6.5 crystal, uj p \ = 1.89 eV @, while for 
5 = 0.09 we have N e g = 0.18 and w p i = 1.68 eV which is 
close to the experimental value. 

Finally, let us compare the relaxation rate T(uj) with 
the optical data for YBCCvs given in Ref. @. In the fre- 
quency region to < 1500 cm -1 and for the temperatures 
T = 0.03i and T = 0.05i, which are close to the exper- 
imental values T = 147 K and T = 244 K, respectively, 
we get the results shown in Fig. [5] As we see, for the op- 
tical properties we also obtain a reasonable quantitative 
agreement of our theory with experiments. 




50 100 150 200 250 300 

T(K) 

FIG. 7: Temperature dependence of the resistivity p(T) at 
(a) S = 0.09 in comparison with the experimental data for 
YBCOe.5 from Ref. Q shown by symbols, and (b) at 8 — 
0.08 (solid) and S = 0.17 (dashed) in comparison with the 
experimental data for LSCO from Ref. [6lJ shown by symbols. 
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1500 



FIG. 8: Relaxation rate T(u) at 5 = 0.09 for T = 0.03t 
(dashed) and T = O.Obt (dotted) in comparison with ex- 
perimental data for YBCO6.5 (Ref. [6J) shown by symbols: 
squares for T = 147 K and triangles for T = 244 K. 



E. Comparison with previous theoretical studies 

Various methods have been used in theoretical studies 
of the optical and dc conductivities in cuprates as dis- 
cussed in Sec. I. Here we compare our results with pre- 
vious studies to clarify what kind of problems the latter 
has and how we have resolved some of them. 

One of the problems is how to explain a linear tem- 
perature dependence of the resistivity in optimally doped 
cuprates in a broad temperature range (see, e.g., Ref. [65 1 
and a discussion in Ref. (s3|). In early studies the 
local-density-functional theory was used in the calcula- 
tion of transport properties of cuprates (for a review 
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see [64j). Calculations of the resistivity within the re- 
laxation rate approximation for electron scattering on 
phonons, (l/rt r ) oc At r T , result in a linear T-dependence 
over a broad temperature range, p cx l/(T tI uj pl ). How- 
ever, the absolute values of the resistivity prove to be 
several times smaller than the experimental ones. This 
discrepancy could be removed by using larger values of 
the transport EPI coupling constant Atr = 1-5 — 2 instead 
of the calculated value A tr = 0.65 — 0.32 for optimally 
doped LSCO and YBCO, or by assuming the theoretical 
plasma frequencies u p i to be smaller by a factor of two 
to three. 

To reconcile a weak transport EPI coupling A tr ^0.5 
and a strong EPI for quasiparticles, \ qp ~ 2, needed 
to explain the high T c in cuprates, a model of a strong 
forward scattering of electrons induced by electron corre- 
lations was proposed 16611 w hich enables to obtain A tr ~ 
X q p/3 (for a review see Sec. III.B). This model was 



used in Ref. 67[ to explain a linear T-dependence of the 
resistivity in a broad region of T where, however, the 
extended van Hove singularities sufficiently close to the 
Fermi energy were assumed and fitting parameters for 
EPI were introduced. 

As shown in Fig. [7] (b), we obtain a linear tempera- 
ture dependence for the resistivity for a nearly optimally 
doped sample in a reasonable agreement with experiment 
on LSCO without using fitting parameters. The values 
of the plasma frequency are also close to experiments. 
In our memory-function theory the transport relaxation 
rate T(u>) and the quasiparticle self-energy are not related 
in a simple way, and therefore we can explain a suffi- 
ciently weak scattering observed in the conductivity and 
a strong superconducting pairing induced by the same 
electron coupling to spin fluctuations (see, e.g., [Hj]). 

To describe the normal state properties, such as the re- 
sistivity and OC, phenomenological spin-fermion models 
have been used. In particular, within a nearly antifcr- 
romagnetic Fermi-liquid model [69| a reasonabl e ag ree- 
ment with experiments on resistivity for YBCO7 |70( and 
on OC for optimally doped and overdoped cuprate com- 
pounds [7l| was obtained by a particular choice of model 
parameters. In Ref. 3 the OC of YBCO*, (x = 6.3, 7) 
within the memory-function method was calculated. To 
obtain an agreement with experiments, several models 
for the spin susceptibility have to be considered. A de- 
tailed qualitative discussion of the OC behavior at vari- 
ous frequencies and temperatures within the spin-fermion 
model was performed in Ref. 42] . Contrary to these phe- 
nomenological approaches, we obtain a fair agreement 
with experiments for various doping and temperatures 
within the microscopic theory for the spin-fluctuation 
susceptibility without fitting coupling parameters. 

There are several studies of the OC within the Hub- 
bard model in the limit of weak correlations, where the 
insulating AF state emerges from a strong AF interac- 
tion as discussed in Sec. I. In particular, in Ref. [73| the 
OC was calculated based on a self-consistent treatment 
of the Hubbard model. It was argued that the charge 



transfer gap observed in the insulating state of cuprates 
is due to the AF LRO which splits the Cu0 2 band into 
two magnetic subbands. Assuming that the AF LRO ex- 
ists at any doping, the experimentally observed increase 
of the charge transfer energy with doping and a simul- 
taneous decrease of the MIR absorption energy were ex- 
plained. In this scenario, the MIR band originates from 
the pseudogap in the electronic spectrum also induced by 
the AF LRO. In our theory based on the consideration 
of the t-J model we cannot study the interesting prob- 
lem of doping dependence of the charge-transfer peak at 
high energy observed in the insulator-to-metal transition 
in cuprates. This problem should be considered within 
the two-subband Hubbard model as, e.g., in Ref. f74j ]. 
since in the t-J model only the lower Hubbard subband 
is explicitly taken into account. 

In Ref. [39j, the MIR absorption in the region of 
lj « t — 2t, being quite strong even without coupling 
to phonons, was related to the interaction of a doped 
hole with spin excitations. In our theory, we can ex- 
plain the MIR absorption by electron interaction with 
spin fluctuations in the decay of charge excitations. We 
cannot relate the MIR absorption to the pseudogap in the 
single-particle electronic spectrum, since the MIR con- 
ductivity shows no notable temperature and doping de- 
pendence which is characteristic for the pseudogap phe- 
nomenon 



75] 



In the limit of strong correlations the Hubbard model 
and t-J model have been used in calculations of dc and 
optical conductivity. As was pointed out in Sec. I, in nu- 
merical studies of finite clusters, due to a finite energy 
resolution, only restricted information on the frequency 
and temperature dependence of the OC can be found. 
The effects of strong correlations have been efficiently 
taken into account within the DMFT method which en- 
ables to reproduce qualitatively the main features of the 
OC: the Drude peak, the MIR region, the charge-transfer 
excitations and the temperature and dop ing dependence 
of the optical spectral weight (see, e.g., [2l|, tZy, l23| and 
references therein). However, in the DMFT the spatial 
correlations, such as short-range AF fluctuations, are not 
taken into account, and therefore the low-energy part of 
the OC caused by charge-boson interaction cannot be 
studied. In the cluster DMFT (28|, due to a finite size of 
clusters, only a qualitative low-frequency behavior of the 
OC can be found. In the present work the complicated 
wave-vector dependence of the dynamical spin suscepti- 
bility and the electron interaction with spin fluctuations 
are fully taken into account without using fitting param- 
eters. This enables to reproduce both the transport re- 
laxation rate T(uj,T) (see Fig. [3]) and the real part of 
the memory function M'(w,T) (see Fig. 0}, yielding the 
optical mass renormalization, in a fair agreement with 
experiments. 

The Allen approximation [40| for the current-current 
correlation function, commonly used in the calculation of 
the OC and the transport relaxation rate 1/t, is based 
on the perturbation theory with respect to (1/lut) <C 1. 
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In our memory-function approach, the optical relaxation 
rate p^|) describes the direct decay of a charge excita- 
tion into an electron-hole pair assisted by the excitation 
of spin fluctuations. In the Allen approach, processes 
of this type appear due to finite life-time effects for the 
electron- hole pair. Therefore, to calculate the optical re- 
laxation rate one has to express the latter in terms of a 
quasiparticle scattering rate which is not a straightfor- 
ward procedure (see, e.g., fiefs. [78|-[8(|). 

The present microscopic theory has in fact some limita- 
tions arising from the t—J model used in the calculations. 
Besides the deficiency of the charge-transfer peak at high 
energy discussed above, the complicated structure of the 
OC found experimentally in the MIR region is missed 
in our theory. This may be due to polaron effects and 
the coupling of magnetic excitations with phonons via 
doped holes as discussed in Refs. 37-39]. To overcome 
these limitations, an extended Hubbard model including 
a strong electron-phonon interaction should be consid- 
ered. 



V. CONCLUSION 

In the present paper we have studied the charge dy- 
namics within a microscopic theory for the optical and dc 
conductivities for the t-J model by taking into account 
electron scattering by spin fluctuations. In our theory, 
based on the memory-function formalism, we calculate 
directly the transport relaxation rate without using the 
Allen perturbation theory. 

Within the proposed theory, we are able to obtain a 
reasonable agreement with experiments on cuprates for 
the relaxation rate, the optical conductivity, and the re- 
sistivity in broad regions of temperatures and doping. 
In particular, in the undcrdopcd region with a strong 
AF SRO, a fair quantitative agreement was found for 
the resistivity, Fig. [7] (a), and for the relaxation rate, 
Fig. [5] This proves the essential role of AF spin fluc- 
tuations in the charge dynamics of cuprates. This con- 
clusion has been corroborated in a number of theoreti- 
cal and experimental studies of OC (see, e.g., @, HI and 
references therein). In the overdoped case, where the 
AF spin fluctuations are suppressed, additional scatter- 
ing mechanisms (e.g., due to electron-phonon interaction 
and impurity scattering) should be invoked to explain ex- 
perimental data. From our results we conclude that spin- 
fluctuations induced by the kinematic interaction should 
give a substantial contribution to the <i-wave pairing in 
cuprates as has been shown recently in Ref. [68j. 
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Appendix A: Calculation of the Memory Function 

To derive Eq. (J6]) for the memory function, we con- 
sider the equations of motion for the relaxation function 
*(t-f) - ((J x (t)\J x (t>))) (see Refs. [H S3, [z3). Dif- 
ferentiating the function subsequently over time t and 
t we obtain a system of equations which in the Fourier 
representation reads, 

w*(w) = Xo + ((F x \J x %, (Al) 
oj ((F X \J X % = -((F X \F X %, (A2) 

where F x = iJ x — [J Xl H] is the force operator. In 
Eq. (|A2[) the relation of the orthogonality (F X ,J X ) = 
(iJ x , J x ) = ([J x , J x \) = was used. Introducing the zero- 
order relaxation function $o( a; ) — Xo/ w we can solve the 
system of equations (IA1[) and (| A2|) in the form 

*(w) = 4 (w)-$o(w)r(w)$ H, (A3) 
with the scattering matrix 

T(w) = (l/xo)((F x \F x )Ul/xo). (A4) 

The memory function M(uS) is defined by Eq. §5§ which 
can be written in the form: 

$(w) = $oM - #o(w) [M M/xo] $(w). (A5) 

From Eqs. (|A3|) and (|A5|) we get a relation between the 
memory function and the scattering matrix: 

T(w) = [M(w)/xo] - [M(w)/xo]$oM2». (A6) 

This equation shows that the memory function is the 
"proper part" of the scattering matrix (|A4I) which has 
no parts connected by a single zero-order relaxation func- 
tion, M(u>) = xoT(uj) plopcl , as given by Eq. ©. 

Appendix B: Mode Coupling Approximation 

Using the spectral representation for the retarded 
Green functions [56[, we write the relaxation rate ©, 
r(w) = M» = Im((F x \F x )Z T °%l(l/ X o), in terms of 
the time-dependent force-force correlation function: 

= ^l_CXp(/^) r°° faiut ( FxFx ( t))P ropcr_ (m) 

2xo^ 

To calculate the force operator F x = [J x ,H], we first 
determine the current J x = —i[P x ,H], where the po- 
larization operator in terms of HOs reads: P x — 
e J2i Rf J2a Xj" 1 ■ Using the commutation relations (fT3|) 
we derive the expression for the current operator 

J x =ieJ2(Ri~ R X j ) Uj XfXf . (B2) 
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The force operator describes electron scattering on spin 
and charge (density) excitations which results from the 
kinematic interaction for the HOs. This can be seen from 
the equation of motion for the electron annihilation op- 
erator: 

^J t Xt(t) = [X?°,H] = -,jX? ct - E UiB^Xf' 
+ (1/2) J ij Xf a '[B jl7r7 , - <W] , (B3) 
where the Bose-like operator Bi aa > is introduced, 

Biaa' — + X^ <T )S cr 'cr + Xf" 5 a i s 

= [1 - (1/2)^ + SflS^a + SJ <W (B4) 

Here the completeness relation for the HOs and the def- 
inition of the number and spin operators (fl4| are used. 
By this type of equations of motion, for the force operator 
we obtain the expression: 

X {X[ \tjlX t a Bjcrcr' — (1/2) Jjl Xj a Blcrcr'} 

-[t tl xf"Bl a , - (1/2) J U xf°Bl a ,]X°°}. (B5) 

Introducing the q-representation for HOs and the inter- 
actions, 

Vn ^ q 

q 

^ = jf E *(<o eiqRi3 . J ^ = jf E J (q) e * qRij - 

q q 

where Rjj = R; — Rj , the force operator (|B5[) takes the 
form, 

Fx = -^EE ^ - (V2)J(q)] 

k,q <j<j' 

^{«- q v-4'Vr s - q ^}, (B7) 

where u x (k) = —dt(k)/dk x is the electron velocity. 
Changing the indexes in the last term, k' = k — q, a ■<-» 



er', and q — >• — q, we obtain the final expression 

F * = ~Jf E Ei ^ k - 1) - (V2) J(q)] 

k,q crcr' 

- u x (k - q) [i(k) - (1/2) J(q)]}^°X^ q B qffCT , 
= - ^ E E ^( k > k - q) *k°*£- q *W ■ (B8) 

k.q crcr' 

In the last equation we introduce the transport vertex 
^(k, k - q) given by Eq. (Q7]). 

We calculate the many-particle time-dependent corre- 
lation functions in Eq. (jBlj) in the mode-coupling approx- 
imation assuming an independent propagation of electron 
and charge-spin excitations. In this approximation, the 
time-dependent correlation functions can be written as a 
product of fermionic and bosonic correlation functions: 

(XZ°X^B^\XZ\(t)X k °(t)Bl a ,(t)) (B9) 
= (X^X^(t))(xt n XZ\(t))(B^B^,(t)) ■ 

Using the definition for the Bose-like operator (|B4[) . for 
the bosonic correlation function we obtain 

(B^Bl a ,(t)) = ({{X™ + X™]6 a , a + X™ <W} 
x{[X°° q (*) + XZ°(fi] K'a + XZ%(t) <W) 
= + X™} {X»\{t) + xz^t)}) 
+(X* a XZl(t)) <W = (l/4)<2V q JV_ q (t)) 5 a - a 

+{S^\SZq{t)) <W + (X™XZl(t)) <W (BIO) 

In the paramagnetic state, for the sum of the spin cor- 
relation functions in Eq. (|B10|) we have: (S' q S'f. q (i)) + 
(SqS-q{t)) = (SqS_q(t)). Finally, using spectral repre- 
sentations for the time-dependent correlation functions 
in Eq. (]B9]), [H 

/OO 
duje- iM f(ij)[-(lM]lm((A\B)) u , 
-OO 

where /(w) is the Fermi function n(uj) for the correlation 
function {X^X^{t)) and the Bose function N(u>) for the 
charge-spin correlation functions, after integration over 
time t in Eq. (|B1[) we obtain the expression (TT5|) for the 
relaxation rate. 



[1] T. Timusk and D. Tanner, in Physical Properties of High 
Temperature Superconductors III, edited by D. M. Gins- 
berg (World Scientific, Singapore, 1992), Chap. 5. 

[2] D.N. Basov and T. Timusk, Rev. Mod. Phys. 77, 721 
(2005). 

[3] D.N. Basov, R. D. Averitt, D. van der Marel, M. Dressel, 



and K. Haule, Rev. Mod. Phys. 83, 471 (2011). 

[4] N.M. Plakida, High- Temperature Cuprate Superconduc- 
tors. (Springer, Heidelberg, 2010) Chap. 5. 

[5] A. A. Kordyuk, V.B. Zabolotnyy, D.V. Evtushinsky, D.S. 
Inosov, T.K. Kim, B. Biichner, and S.V. Borisenko, Eur. 
Phys. J. Special Topics 188, 153 (2010). 



12 



J. Hwang, J. Yang, T. Timusk, S.G. Sharapov, J. P. Car- 
botte, D.A. Bonn, R. Liang, and W.N. Hardy, Phys. Rev. [31 
B 73, 014508 (2006). [32 
J. Yang, J. Hwang, E. Schachinger, J. P. Carbotte, 
R.P.S.M. Lobo, D. Colson, A. Forget, and T. Timusk, [33 
Phys. Rev. Lett. 102, 027003 (2009). 

E. van Heumen, E. Muhlethaler, A. B. Kuzmenko, H. [34 
Eisaki, W. Meevasana, M. Greven, and D. van der Marel, 
Phys. Rev. B 79, 184512 (2009). [35 
Ph. Bourges, in: The Gap Symmetry and Fluctuations 
in High Temperature Superconductors, ed. by J. Bok, [36 
G. Deutscher, D. Pavuna, and S. A. Wolf (Plenum Press, [37 
1998), p. 349-371. 

M. Le Tacon, G. Ghiringhelli, J. Chaloupka, M. Moretti [38 
Sala, V. Hinkov, M. W. Haverkort, M. Minola, M. Bakr, 
K. J. Zhou, S. Blanco-Canosa, C. Monney, Y. T. Song, 
G. L. Sun, C. T. Lin, G. M. De Luca, M. Salluzzo, G. [39 
Khaliullin, T. Schmitt, L. Braicovich, and B. Keimer, 
Nature Phys. 7, 725 (2011). [40; 
E.G. Maksimov, M.L. Kulic, and O.V. Dolgov, Adv. in [41 
Cond. Mat. Physics, DOI: 10.1155/2010/423725 (2010). 
S. Uchida, T. Ido, H. Takagi, T. Arima, Y. Tokura, and [42 
S. Tajima, Phys. Rev. B 43, 7942 (1991). 

S. Uchida, K. Tamasaku, and S. Tajima, Phys. Rev. B [43 
53, 14558 (1996). 

Y. Onose, Y. Taguchi, K. Ishizaka, and Y. Tokura, Phys. 
Rev. B 69, 024504 (2004). 

Y.S. Lee, K. Segawa, Z.Q. Li, W.J. Padilla, M. Dumm, [44 
S.V. Dordevic, C.C. Homes, Y. Ando, and D.N. Basov, [45 
Phys. Rev. B 72, 054529 (2005). 

W.J. Padilla, Y.S. Lee, M. Dumm, G. Blumberg, S. Ono, [46 
K. Segawa, S. Komiya, Y. Ando, and D.N Basov, Phys. [47 
Rev. B 72, 060511(R) (2005). [48 
S. Lupi, D. Nicoletti, O. Limaj, L. Baldassarre, M. Or- 
tolani, S. Ono, Yoichi Ando, and P. Calvani, Phys. Rev. [49 
Lett. 102, 206409 (2009). 

Y. Iye, in Properties of High Temperature Superconduc- [50 
tors, edited by D. M. Ginsberg (World Scientific, Singa- [51 
pore, 1992), Vol. 3, pp. 285 - 361. [52 
J. Zaanen, G.A. Sawatzky, and J.W. Allen, 
Phys. Rev. Lett. 55, 418 (1985). [53 
P.A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. [54 
78, 17 (2006). 

A. Comanac, L. De' Medici, M. Capone, and A.J. Millis, 
Nature Physics 4, 287 (2008). [55 
A. Georges, B. G. Kotliar, W. Krauth, and M. J. Rozen- [56 
berg, Rev. Mod. Phys. 68, 13 (1996). [57 
C. Weber, K. Haule, and G. Kotliar, Phys. Rev. B 78, 
134519 (2008). [58 
C. Weber, K. Haule, and G. Kotliar, Nature Physics 6, 
574578 (2010). [59 
G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 
O. Parcollet, and C. A. Marianetti, Rev. Mod. Phys. 78, 
865 (2006). [60 
N.M. Plakida and V.S. Oudovenko, JETP 104, 230 
(2007). [61 
K. Haule and G. Kotliar, Phys. Rev. B 76, 104509 (2007). 
K. Haule and G. Kotliar, Europhys Lett. 77, 27007 
(2007). 

[29] F. Carbone, A.B. Kuzmenko, H.J.A. Molegraaf, E. van [62 
Heumen, V. Lukovac, F. Marsiglio, D. van der Marel, K. 
Haule, G. Kotliar, H. Berger, S. Courjault, P.H. Kes, and [63 
M. Li, Phys. Rev. B 74, 064510 (2006). 

[30] N. Bontemps, R.P.S.M. Lobo, A.F. Santander-Syro, and [64 



A. Zimmers, Ann. Phys. 321, 1547 (2006). 

E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

W. Stephan and P. Horsch, Int. J. Mod. Phys. B 6, 589 

(1992). 

T. Tohyama, P. Horsch, and S. Maekawa, Phys. Rev. 
Lett. 74, 980 (1995). 

R. Eder, Y. Ohta, and S. Maekawa, Phys. Rev. Lett. 74, 
5124 (1995). 

J. Jaklic and P. Prelovsek, Advances in Phys. 49, 1 
(2000). 

M. Vojta and KW. Becker, Eur. Phys. J. B 3, 427 (1998). 

E. Cappelluti, S. Ciuchi, and S. Fratini, Phys. Rev. B 76, 
125111 (2007). 

A.S. Mishchenko, N. Nagaosa, Z.-X. Shen, G. De Filippis, 
V. Cataudella, T.P. Devereaux, C. Bernhard, KW. Kim, 
and J. Zaanen, Phys. Rev. Lett. 100, 166401 (2008). 
L. Vidmar, J. Bonca, and S. Maekawa, Phys. Rev. B 79, 
125120 (2009). 

P.B. Allen, Phys. Rev. B 3, 305 (1971). 

S.V. Shulga, O.V. Dolgov, and E.G. Maksimov, Physica 

C (Amsterdam) 178, 266 (1991). 

Ar. Abanov, A.V. Chubukov, and J. Schmalian, Ad- 
vances in Phys. 52, 119 (2003). 

A.V. Chubukov, D. Pines, and J. Schmalian, in: The 
Physics of Conventional and Unconvencional Supercon- 
ductors, ed. by K.H. Bennemann and J.B. Ketterson, 
(Springer Verlag, Heidelberg-Berlin, 2004) Vol. I, p. 495. 
H. Mori, Prog. Theor. Phys. 34, 399 (1965). 
W. G6tze and P. Woffle, J. Low Temp. Phys. 5, 575 
(1971); Phys. Rev. B 6, 1226 (1972). 
D. Ihle and N.M. Plakida, Z. Phys. B 96, 159 (1994). 
N.M. Plakida, Z. Phys B 103, 383 (1997). 
A. A. Vladimirov, D. Ihle, and N.M. Plakida, Phys. Rev. 
B 80, 104425 (2009). 

A. A. Vladimirov, D. Ihle, and N.M. Plakida, Phys. Rev. 
B 83, 024411 (2011). 

R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 
P.W. Anderson, Science 235, 1196 (1987). 

F. C. Zhang and T.M. Rice, Phys. Rev. B 37, 3759 
(1988). 

J. Hubbard, Proc. Roy. Soc. A 285, 542 (1965). 
Yu.A. Izyumov and Yu.N. Scryabin, Statistical Mechan- 
ics of Magnetically Ordered Systems, (Consultant Bu- 
reau, New York, 1989) pp.74 - 108. 

F. Dyson, Phys. Rev. 102, 1217 (1956). 

D.N. Zubarev, Sov. Phys. Uspekhi 3, 320 (1960). 

G. Jackeli and N.M. Plakida, Phys. Rev. B 60, 5266 
(1999). 

N.M. Plakida and V.S. Oudovenko, Phys. Rev. B 59, 
11949 (1999). 

A. A. Kordyuk, S.V. Borisenko, A. Koitzsch, J. Fink, 
M. Knupfer, and H. Berger, Phys. Rev. B 71, 214513 

(2005) . 

P. Prelovsek and A. Ramsak, Phys. Rev. B 65, 174529 
(2002). 

G.S. Boebinger, Y. Ando, A. Passner, T. Kimura, 
M. Okuya, J. Shimoyama, K. Kishio, K. Tamasaku, 
N. Ichikawa, and S. Uchida, Phys. Rev. Lett. 77, 5417 
(1996). 

Y. Ando, S. Komiya, K. Segawa, S. Ono, and Y. Kurita, 
Phys. Rev. Lett. 93, 267001 (2004). 

Y.C. Ma and N.L. Wang, Phys. Rev. B 73, 144503 

(2006) . 

W.E. Pickett, Rev. Mod. Phys. 61, 433 (1989). 



13 



[65] T. A. Friedmann, M.W. Rabin, J. Giapintzakis, J. P. 

Rice, and D.M. Ginsberg, Phys. Rev. B 42, 6217 (1990). 
[66] R. Zeyher and M. L. Kulic, Phys. Rev. B 53, 2850 (1996). 
[67] G. Varelogiannis and E. N. Economou, Europhys. Lett. 

42, 313 (1998). 

[68] N.M. Plakida, e-arXiv: cond-mat/1110.6715 (to be pub- 
lished in: J. Supercond. Nov. Magn. (2012)). 

[69] A. J. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 
167 (1990). 

[70] P. Monthoux and D. Pines, Phys. Rev. B 49, 4261 (1994). 
[71] B. P. Stojkovic and D. Pines, Phys. Rev. B 56, 11931 
(1997). 

[72] B. Arfi, Phys. Rev. B 45, 2352 (1992). 
[73] T. Das, R.S. Markiewicz, and A. Bansil, Phys. Rev. B 
81, 174504 (2010). 



[74] N.M. Plakida, J. Phys. Soc. Jpn. 65, 3964 (1996). 

[75] S. Hiifner, M.A. Hossain, A. Damascelli, and G.A. 

Sawatzky, Rep. Prog. Phys. 71, 062501 (2008). 
[76] M. Jarrell, J. K. Freericks, and Th. Pruschke, Phys. Rev. 

B 51, 11704 (1995). 
[77] A. Toschi, M. Capone, M. Ortolani, P. Calvani, S. Lupi, 

and C. Castellani, Phys. Rev. Lett. 95 097002 (2005). 
[78] J. P. Carbotte, E. Schachinger, and J. Hwang, 

Phys. Rev. B 71, 054506 (2005). 
[79] E. Schachinger, D. Neuber, and J. P. Carbotte, Phys. 

Rev. B 73, 184507 (2006). 
[80] S.G. Sharapov and J. P. Carbotte, Phys. Rev. B 72, 

134506 (2005). 



